长江下游地区50-80天降水低频分量次季节预测模型实验数据集

杨秋明

江苏省气象科学研究所，南京 210009

摘 要：运用1979-2014年长江下游逐日降水量数据（25个站点的平均值）和NCEP/NCAR全球850 hPa纬向风格点再分析数据（2.5°×2.5°），经过经验正交函数（EOF）和奇异谱分析（SSA），得到用于构建扩展复数自回归模型（ECAR）的次季节预测实验数据集。对2001-2014年长江下游降水低频分量进行次季节逐日变化回报试验的结果表明，50-80 d时间尺度的长江下游低频降水分量的预测时效可达52 d左右，预报能力明显优于自回归模型（AR），而且6-8月的预报技巧最高。该数据集包括1979年1月1日-2014年12月31日全球850 hPa 50-80 d低频纬向风场的7个低频主成分pc1，pc2，…，pc7和长江下游地区逐日低频降水（the daily low frequency rainfall in the lower reaches of the Yangtze River Valley，lcjr）的时间序列（序列长度是13,149 d）。数据集存储在Excel文件，数据量1.19 MB。

关键词：长江下游；降水；50-80天低频分量；次季节预测

DOI: https://doi.org/10.3974/geodp.2018.02.07

1 前言

近10年来，气候观测获得的全球数据迅速增加，这些具有高度数据相关性和多重数据属性的科学大数据反映和表征着复杂的自然现象与关系[1–2]。从大量数据中提取出部分有效数据，可以获得较全面的大气季节内振荡（ISO）变化信息，为极端天气事件15-60 d次季节预报[3–5]提供了较好的发展基础。杨秋明（2016）基于全球低频环流和长江下游低频降水数据，建立了简化的时变扩展复数自回归预测模型（ECAR）[6]，明显延长了预测时效，可以较好地预测未来50-60 d内与长江下游地区强降水显著相关的50-80 d降水低频分量的变化[7]。这些次季节预测数据来源于NCEP/NCAR全球850 hPa纬向风格点再分析数据（2.5°×2.5°）和1979-2014年长江下游逐日降水量数据（30°30′N-32°0′N，118°0′E-122°30′E；25个站点的平均值），通过带通滤波、EOF和SSA等分析，得到1979年1月1日-2014年12月31日全球850 hPa 50-80 d低频纬向风场的7个低频主成分pc1，pc2，…，pc7和长江下游地区逐日低频降水lcjr的时间序列，构成长江下游地区50-80 d降水低频分量次季节预测数据集。

2 数据集元数据简介

长江下游地区50-80 d降水低频分量次季节预测数据集[8]的名称、作者、地理区域、数据年代、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。

表1 长江下游地区50-80 d降水低频分量次季节预测数据集元数据简表
条目描述
数据集名称长江下游地区50-80 d降水低频分量次季节预测数据集
数据集短名ForecastLowFreqRainfallLYRV
作者信息杨秋明 G-9579-2018，江苏省气象科学研究所，yqm0305@263.net
地理区域中国长江下游地区，30°30′N-32°0′N，118°0′E-122°30′E
数据年代1979-2014年
数据格式.xls
数据量1.19 MB
数据集组成1979-2014年逐日全球850 hPa纬向风低频主成分和长江下游50-80天低频降水
基金项目国家自然科学基金项目（41175082）
出版与共享服务平台全球变化科学研究数据出版系统 http://www.geodoi.ac.cn
地址北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所
数据共享政策全球变化科学研究数据出版系统的"数据"包括元数据（中英文）、实体数据（中英文）和通过《全球变化数据学报》（中英文）发表的数据论文。其共享政策如下：（1）"数据"以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用"数据"需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）"数据"的用户需要与《全球变化数据学报》（中英文）编辑部签署书面协议，获得许可；（4）摘取"数据"中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[9]

3 数据研发方法

850 hPa纬向风场选用NCEP/NCAR逐日再分析全球风场资料（2.5°×2.5°格点）[10]，逐日降水量采用长江下游地区（30°30′N-32°0′N，118°0′E-122°30′E）25个站的平均值，以上资料时间是1979年1月1日-2014年12月31日。

3.1 算法原理

对全球低频纬向风进行主成分分析（PCA），得到前7个主要空间模态（时间尺度是50-80 d，资料时间是1979-2000年，序列长度是 8,036 d，由Butterworth滤波得到全球850 hPa低频纬向风场），解释方差分别是23.5%、4.3%、4.1%、3.5%、3.0%、2.9%和2.5%[7]。用奇异谱分析（Singular Spectrum Analysis，SSA）[11]，对长江下游逐日降水原始序列和全球纬向风主要模态的时间系数进行低通滤波，重建对应于主要的月以上时间尺度次季节振荡信号的分量序列，得到观测的长江下游降水低频分量序列和全球850 hPa逐日纬向低频风场空间模态对应的低频主成分。

3.2 技术路线

将2001-2014年逐日观测的全球纬向风投影到这7个低频空间模态（由1979-2000年资料计算，序列长度是8,036 d），得到前7个主成分pc1-pc7的观测值（包含逐日高频扰动）；然后将这7个观测的pc和同期的长江下游逐日降水投影到各自经SSA后的对应于50-80 d振荡的T-EOF上得到低频主成分pc1，pc2，…，pc7和长江下游逐日降水的50-80 d重构分量lcjr（序列长度是5,113 d）。

4 数据结果与验证

4.1 数据集组成

该数据集包括1979年1月1日-2014年12月31日全球850 hPa 50-80 d低频纬向风场的7个低频主成分pc1，pc2，…，pc7和长江下游地区逐日低频降水lcjr的时间序列（序列长度是13,149 d）。数据集存储在Excel文件，数据量1.19 MB。

4.2 数据结果

数据结果包括：

（1）1979年1月1日-2014年12月31日全球850 hPa 50-80 d低频纬向风场的7个低频主成分pc1，pc2，…，pc7；

（2）1979年1月1日-2014年12月31日长江下游地区逐日50-80 d低频降水lcjr的时间序列。

4.3 数据结果验证

2013-2014年长江下游50-80d低频降水与全球850 hPa低频纬向风主成分的200d滑动相关

图1 2013-2014年长江下游50-80 d低频降水lcjr与全球850 hPa低频纬向风主成分pc1-pc7(a-g)的200 d滑动相关

下载

采用滑动建模的方法，使得长江下游低频降水lcjr与pc1，pc2，...，pc7的各个时段的相关中尽可能保持在较显著状态，逐步适应观测数据中的主要时滞相关变化，从而提高ECAR[6]的预测能力。文中的实时预测用限定记忆法，保持子序列N0不变，由扩展数据阵F = (pc1 pc2 ... pc7 lcjr)，构建含时变系数的ECAR对扩展数据阵中的各个数据分量滑动进行独立样本预测试验。图1给出了2013年4月-2014年8月期间长江下游50-80 d降水低频分量与全球低频纬向风主要模态对应的时间系数pc1，...，pc7之间200 d滑动相关的时间演变。结果分析表明，滑动相关系数均具有显著的季节和年际变化，强正负相关时段具有季节性和不规则的时间变化。仅与pc3（pc7）的滑动相关表现为较稳定的正（负）相关（图1c，1g），其余的均呈现相关符号改变的不稳定，尤其是与pc4、pc5、pc6的滑动相关（图1d，1e，1f）分别在2014年4月、2013年7月和2013年10月附近发生显著正负（或负正）相关的突变，与全球环流50-80 d低频分量与东亚季风降水之间的非线性相互作用及其时间变化有关。因此，这种建模方案可以适应低维空间中各个低频分量之间的联系随时间的改变，反映气候系统分量之间相关的季节和年际变化。

5 讨论和总结

气象历史和实时大数据可以代表全局，揭示出全局特有的宏观信息，这些综合信息是过去较小规模数据难以挖掘的。通过数据分解、变换，直接从多变量、超高维、超长序列、高度耦合相关和高度时空属性的观测资料中提取一定时间尺度的、与一定区域（如长江流域等）的极端天气气候事件关系密切的各种主要低频振荡型，研究其不同的传播特征和强度的季节和年际变化规律，由动态数据来驱动复杂低频变化过程与系统的构建，建立简化的次季节预测模型，可以显著延长尺度低频分量和特定区域中对应的极端天气气候事件预报时效。所以，我们研制了全球850 hPa 50-80 d低频纬向风场的7个低频主成分pc1， pc2，…，pc7和长江下游地区逐日低频降水lcjr数据集，用于建立含时变系数的ECAR预测模型。需要注意的是50-80 d低频振荡强度显示出不规则的时间变化，因此后续工作需要将时间序列延长，及时更新。

参考文献